arXiv: 1504.023 lOvl [physics.bio-ph] 9 Apr 2015 


Zipping and Entanglement in Flagellar Bundle of E. Coli: Role of Motile Cell Body 


Tapan Chandra Adhyapak^’* and Holger Stark^A 

^ Institut fur Theoretische Physik, Technische Universitdt Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany 

(Dated: April 10, 2015) 

The course of a peritrichous bacterium such as E. coli crucially depends on the level of synchro¬ 
nization and self-organization of several rotating flagella. However, the rotation of each flagellum 
generates counter body movements which in turn affect the flagellar dynamics. Using a detailed 
numerical model of an E. coli, we demonstrate that flagellar entanglement, besides fluid flow relative 
to the moving body, dramatically changes the dynamics of flagella from that compared to anchored 
flagella. In particular, bundle formation occurs through a zipping motion in a remarkably rapid 
time, affected little by initial flagellar orientation. A simplified analytical model supports our obser¬ 
vations. Finally, we illustrate how entanglement, hydrodynamic interactions, and body movement 
contribute to zipping and bundling. 


Understanding self-propulsion of microorganisms pose 
utmost challenges involving rich and complex physics [1- 
7]. Bacteria are among the simplest and widely stud¬ 
ied of such systems [8-13]. Yet, only recently we are 
able to explore in full detail the underlying physics in¬ 
volved [14-26]; much of this, however, is still to be ap¬ 
prehended. Propulsion of peritrichous bacteria such as 
E. coli is generated by the rotation of a bundle of several 
helical propellers, called flagella. Flagella are passive fil¬ 
aments rotated at one end by rotary motors embedded 
in the cell wall [27]. The level of synchronization and 
self-organization of rotating flagella crucially decides the 
swimming course of the cell body, to which they are at¬ 
tached, making it either propel or tumble. However, the 
dependency is not one-sided: the rotary motors that ro¬ 
tate each flagellum also produce systematic body move¬ 
ments, which in turn affect the flagellar dynamics. While 
the rotating cell body drags the proximal ends of flagella 
with it, the distal ends cannot keep up due to friction 
with the surrounding fluid. In understanding flagellar 
synchronization and bundling dynamics, focus has so far 
been given primarily on hydrodynamic interactions and 
elastic properties of flagella [19, 21, 25, 26, 28, 29]. Al¬ 
though, body movement is speculated to play an impor¬ 
tant role too [19, 30], knowledge of its detailed impact is 
still lacking. 

The cell body moves in response to the forces and 
torques acting on it [31]. It translates due to the thrust 
force generated by the rotating bundle of flagella in the 
surrounding fluid medium at low Reynolds number. It 
also has to rotate since the torque on the body has to 
balance all motor torques acting on the flagella. Ref. [30] 
argued that the sole effect of body rotation on a flagellum 
is to simply wrap it around the cell-body axis and thereby 
enhance bundling. However, real flagella are helical, dur¬ 
ing tumbling they are more or less arbitrarily oriented 
[8], and they cannot simply pass through each other. We 
will demonstrate through our simulations that body ro¬ 
tation in such situations leads to entanglement, where 
portions of different flagella obstruct each other’s free 
course. While it is well established that hydrodynamic 


interactions are sufficient to synchronize [26] and bundle 
[18, 19, 21] anchored flagella, the role of entanglement 
for flagellar dynamics is not known so far. In particular, 
an understanding how cell body motion influences flagel¬ 
lar dynamics is incomplete without considering flagellar 
entanglement. 

In this paper we explore how body movements influ¬ 
ences the dynamics of flagella with the help of a realis¬ 
tic numerical model of an E. coli that includes detailed 
flagellar elasticity, hydrodynamic and steric interactions 
among flagella, and a motile cell body. We demonstrate 
that body movements dramatically change flagellar be¬ 
havior leading to profound impacts on the overall dynam¬ 
ics of the cell. In particular, bundle formation happens 
through flagellar synchronization and a ‘zipping’ motion 
on an experimental time scale [8] , which is much smaller 
than for anchored flagella [19, 21] and which is approxi¬ 
mately independent of the initial orientations of unbun¬ 
dled flagella. We also analyze the relative importance of 
body movements and flagellar interactions for synchro¬ 
nization and bundling. Our work therefore is a major 
step towards understanding the propulsion of a peritric¬ 
hous bacterium close to its real conditions [8, 9]. 

We first summarize our approach to describe the dy¬ 
namics of the cell body with multiple flagella in an 
unbounded fluid of viscosity ij. We treat each flagel¬ 
lum as a slender body with centerline r(s) parametrized 
by the arc length s. By affixing the orthonormal tri¬ 
pod {ei(s), 62(5), 63(5)} at each point on the centerline, 
where 03 is the local tangent and ei and 82 are unit vec¬ 
tors along the principal axes of the flagellar cross section, 
one can fully characterize the bent and twisted flagellum. 
Dynamics of the flagellum now is governed by Langevin 
equations for r(s) and the twist angle <j){s) about the 
centerline [17]: 

dt^ = Pt (Fei + Fg -1- Fth) + vjj, (1) 

dt(l) = Pr{Tel + Tth), (2) 

Here, we separate the velocity contribution Vh due to 
hydrodynamic interactions from local terms and denote 
local forces and torques by F’s and T’s, respectively. 


2 



FIG. 1. Schematics of the cell body and part of one flag¬ 
ellum. For clarity, length h between discretization points 
Ti is out of proportion and only the first material tripod 
{ei(l), 62 ( 1 ), 63 ( 1 )} is shown. The cell body is centered at vi, 
and oriented along the vector 653 . The flagellum is attached 
at To off-centered by a distance d along ef, 3 . The motor torque 
drives a tripod at ro. 

Self-mobilities = 63 0 es/qy -|- (I — 03 (g) e 3 )/ 7 _L and 
= l/ 7 it are expressed in terms of friction coefficients 
per unit length. For the flagellum of an E. coli they 
are yy = 1.6 x 10“^pNs//im^, 7 _l = 2.8 x 10“^pNs/^m^ 
and 7 fl = 1.26 x 10“®pNs [20]. Thermal forces Fth and 
torques Tth are shown for completeness. Although these 
are predicted to play an important role during rotation- 
induced polymorphic transformations of a flagellum [ 20 ] , 
they are negligible in our present study [ 21 ] and are ig¬ 
nored. 

Elastic forces and torques Fei = —SE/6r and Tei = 
—SEjScj) are derived, respectively, from the total elas¬ 
tic free energy [r(s), (/)(s)] of the flagellum, the form 
of which is obtained as follows. The rotational strain 
vector El moves the material tripod along the flagellum: 
dsBi, = El X V = 1,2,3. Therefore, its components 
completely characterize the instantaneous flagellar con¬ 
formation [20] . A small deformation dEl = — $7° from 

the normal helical ground state = {0.0,1.3,—2.1} 
[17, 32] needs the Kirchhoff elastic free energy density 
[33], /k(f2) = (A/2) [(dfli)2 + {dQ2?] + {C/2){dn^f. 
For E. coli, we choose an isotropic bending rigidityA = 
3.5pN/j,m^ (assuming a circular flagellar cross section) 
and the twist rigidity C = A [9]. Integrating over 
the length of the flagellum, we obtain [r(s), (/(s)] = 
/ ds (/k + fSt), where we include a stretching free energy 
density fst = K{dsr)'^/2 with K = lO^pN [17]. 

To proceed, we discretize Eqs. (1) and (2) by consider¬ 
ing discrete positions = r(si) along each flagellum and 
by assigning {ei(l), 02 ( 1 ), 63 ( 1 )} to the straight segment 
of length h between Vi_i and (see Fig. 1). Excluded- 
volume interactions among flagella are enforced by the 
steric force Fs(ri) = - F{(/i — hj)/h. Here, the sum¬ 

mation runs over all overlaps occurring within [r^-i, r^+i] 
of a given flagellum and F{ is the steric force at a distance 
hj < h from r^, appropriately decomposed to act on the 
discrete points (for details see the supplemental mate¬ 
rial [34]). F{ derives from the Lennard-Jones potential 
ULjirm) = (Foa/G) [(cr/r„)^^ - {a/rmf] Q{2^/^a - r^). 
We truncate it at the mimimum using the heaviside step 


function 0 (a;), where Vm is the minimal distance be¬ 
tween the approaching two flagellar centerlines and Fq 
is the strength of the steric force at Vm = cr- We choose 
Fq = 0.8 pN and adjust tr = 4a, with a the cross-sectional 
radius of the flagellar filament, ensuring numerical sta¬ 
bility during entanglement in all situations. 

Finally, to include hydrodynamic interactions between 
the flagella, we treat each discrete point as a sphere of 
radius a and set Vh(ri) = Here, is the 

Rotne-Prager mobility matrix [35] for spheres at and 
Tj, F(rj) is the local force at r^, and the summation runs 
over all points of both flagella. We neglect subleading 
effects from hydrodynamic interactions due to rotation 
of the spheres. Furthermore, neglecting hydrodynamic 
lubrication for close flagella is justified for thin filaments 
and the presence of asperities in real flagella [36]. 

We model the cell body by a spherocylinder of length 
Lh = 2.5/rm and width d), = 0.8/rm [27] (see Fig. 1). 
Point ro of each flagellum is fixed on the body surface 
and a motor torque = T^m drives the flagellum 
by rotating the motor tripod { 01 ( 0 ), 02 ( 0 ), 03 ( 0 ) = m) 
at ro. This tripod couples to the main part of the flagel¬ 
lum through the Kirchhoff elastic free energy density /k, 
where we set A —>• 0 and C —^ 3(7. Thus, the driving 
torque is transferred to the flagellum through a ‘hook’ 
that acts as an universal joint with low bending and high 
twist rigidities [37] allowing the first flagellar segment 
along 03 ( 1 ) to be at any angle to m. In response, the 
body moves and rotates with velocities and 

u>b = M 6 (T& + Tm), respectively. Here, F^ and T;, are the 
respective force and torque (relative to the body center) 
resulting from the forces Fei -t- Fg that act on the flagellar 
anchoring points. For the mobilities /x} and /x} we use 
the analytically available values for a prolate spheroid of 
aspect ratio Lb/db [38]. The angle (j)m between m and 
0 t ,3 is expected to differ from 90° because of, for exam¬ 
ple, a locally curved body surface. We adjust (j)m = 55° 
to obtain a ratio for the bundle-to-body rotation within 
the experimentally observed range [9]. Furthermore, we 
employ the same potential ULji'i’m) to describe the ex¬ 
cluded volume interaction between the body and flagella, 
where Vm now is the minimal distance between the body 
surface and any ri on a flagellum. 

We now turn to our observations. Fig. 2 shows typi¬ 
cal snapshots of the bacterium moving towards the left 
with the flagella in their normal left-handed helical form. 
The snapshots were obtained at regular intervals from 
a particular simulation run. Quantitative details of the 
corresponding flagellar dynamics are presented in Fig. 3. 
At time t = 0 ms, flagella start with an angle between 
their axes. As time progresses, they rotate counterclock¬ 
wise about their axes (as viewed from behind the cell) 
when driven by a positive motor torque Tm = 3.4 pN/rm 
[9, 17]. Simultaneously, the cell body performs a counter¬ 
balancing clockwise rotation and also translates because 
of the thrust force generated by the flagella. The resul- 





3 



FIG. 2. As the bacterium moves to the left, the flagella form 
a bundle of length G similar to a closing zipper. Snapshots 
at various times t (in ms) for <1? = 90° are shown. Rotation 
directions of cell body and flagella are depicted. 


tant flagellar evolution is complex. It involves entangle¬ 
ment and large bending of flagellar axes. To quantify 
synchronization of flagellar rotation in such situations, 
we need to compare the respective tripod vectors e^{i) 
and e^{i) from the two flagella at the same flagellar posi¬ 
tion i. Therefore, we introduce the effective phase differ¬ 
ence 9{i) = cos“^ [ei(*) • /tt with e*(z) = 7?.[ei(i)], 
where TZ rotates the tangent vectors onto each other 
[ 63 ( 1 ) = 63 ( 1 )] about the axis 63 ( 1 ) x 63 ( 1 ). Start¬ 
ing from a non-zero initial value, the contour average 
(6) quickly drops towards zero and the flagella reach a 
nearly synchronized rotational state after about 10 ms 
[see Fig. 3(a)]. This is also reflected in the snapshot at 
10 ms shown in Fig. 2, where the phases of both flag¬ 
ella clearly match. The initial regime remains unaffected 
when changing the initial value of {0). It is completely 
determined by hydrodynamics since steric forces are zero 
as documented by the inset of Fig. 3(a). The subsequent, 
almost linear decrease of (9) towards full synchronization 
coincides with bundle formation, which we discuss now. 

While cell-body movements only quantitatively change 
synchronization dynamics, they dramatically influence 
bundling dynamics compared to the case of anchored 
flagella [18, 19, 21 ]. At t ^ 10ms proximal portions of 
the flagella start to wrap around the body axis [see Fig. 
2 ] , after synchronization has already proceeded consider¬ 
ably. This marks the beginning of bundling near the cell 
body, while rest of the flagella are still apart. With time 
the front of the bundled portion advances away from the 
body, gradually drawing remaining loose flagellar por¬ 
tions into the bundle. The growth of the bundled portion 
thus resembles a ‘zipping’ motion, where the ‘zip’ starts 
at a point near the body and continues till whole of the 
flagella have joined the bundle. 

To quantify these findings, we define the bundled por¬ 
tion as the part of a flagellum for which all its points are 
at distances < 2R from the other flagellum. We choose 
R = 0.22/rm to be the equilibrium radius of the helical 
flagellum. The bundle length A is then measured from 
the cell body to the bundle front normalized by the axial 
length Lax of a flagellum. As seen from Fig. 3(b), ‘zip¬ 
ping’ takes place within approximately 5-20 ms. During 
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FIG. 3. (a) Gontour-averaged phase difference {9} plotted 

versus time t for various values of the opening angle <I>. (b) 
Bundle length L plotted against t for the same 4>'s obtained 
from simulations (full lines) and an analytical model (dashed 
lines). The analytical model fits best for R = 0.30pm, when 
Vc = 22.5pm/s, LOcl'l'n = 11.4 Hz, A/ = 2.5pm, Vf = 83.3 
Hz. (c) Body motion (vc, Wc) and flagellar phase velocity 
Ve. determine the speed of the bundle front P in the analytic 
model. Insets; Steric force density /c (in pN/pm) correlated 
with time evolutions of (fi) in (a) and R in (b), for 4> = 44° 
and 90°. 


this period, after an initial sharp increase A grows al¬ 
most linearly with time. Furthermore, for the studied 
range of opening angles $ the bundling time only varies 
by about 10 ms. Therefore, it is always much smaller than 
the total tumble time of 150 - 450ms [ 8 ]. So, bundling 
during a tumble event is remarkably independent of the 
extent to which a flagellum is thrown out of a bundle. 
This gives a crucial insight into the locomotion of a bac¬ 
terium since bundling without supporting body move¬ 
ments takes much longer. 

Flagellar entanglement is observed to play an impor¬ 
tant role in all the findings mentioned above. During 
zipping, the steric force density fc = [F^l/L, where L 
is the length of a flagellum, is found to build up until 
bundling completes, after which fc declines rapidly [in¬ 
set, Fig. 3(b)]. Due to viscous drag and flagellar flexibil¬ 
ity, parts of the flagella not in the bundle hardly follow 
the cell-body rotation (see movie Ml in the Supplemen¬ 
tal material [34]). As a result, proximal portions of the 
flagella not only start to wrap around the body axis but 
also get entangled [snapshot at f = 10 ms, Fig. 2]. This 
drastically enhances flagellar bundling: further rotation 
of proximal ends is possible only when the entangled front 
proceeds away from the body, gradually bringing the rest 
of the flagella quickly into the bundle. A simplified ana- 



























4 


lytical model discussed below further supports these ob¬ 
servations. A correlation between fc and the slow linear 
decline of {9) mentioned earlier is clear from the inset of 
Fig. 3(a). While flagella in the bundle are synchronized, 
the local phase differences 9{i) of their free ends fluctuate 
strongly until entanglement forces them into the bundle. 
This is observed to be responsible for the delayed slow 
decline of {9). 

A simplified model for the zipping dynamics takes into 
account flagellar entanglement at the bundle front and 
fluid flow in the body-fixed reference frame. The latter 
occurs with respective translational and rotational veloci¬ 
ties Vc = — Vf, and ujc = —Wb [see Fig. 3(c)]. Accordingly, 
the length of the bundle grows with speed dlr/dt and the 
angle ip = 90° — $/2 varies in time according to 


dir ^cR 

- = Ue H- 

dt cos Ip 


, dip Vc cos Ip 

and -j- = -f - y 

^ax 


( 3 ) 


Like in a zipper the free portions of the rotating flag¬ 
ella are dragged into the bundle front P with the helical 
phase velocity Vg. They are perfectly fit into the bundle; 
hence Ir grows with the same speed Vg = where 

A/ is the helical pitch and Vf the frequency of rotation. 
Second, the surrounding fluid wraps the free flagellar por¬ 
tions with an angular velocity uig onto the bundle cylin¬ 
der. So, flagella are dragged with an additional speed 
uJcR/ cosip into the bundle front, where R is an effective 
bundle radius. Finally, the translational flow rotates the 
flagellar tip with speed Vg cos ip about the bundle front 
P, which gives the angular velocity dip/dt. Eqs. (3) are 
solved for lr(t) using parameter values measured from 
our simulations with R and lr{t = 0) adjusted for the 
best fit. The results plotted as dashed lines in Fig. 3(b) 
agree well with our simulation results, strengthening the 
interpretation of our observations discussed above. 

To obtain further insight into flagellar dynamics, we 
make a comparative study to judge the relevance of cell- 
body movements and various flagellar interactions. The 
results are presented in Fig. 4. We consider the stan¬ 
dard simulation run performed at d* = 23°. Its results 
[curves (i) in Fig. 4] are compared with those obtained 
from new simulations where either body movements or 
flagellar interactions are ignored. The impact of body 
rotation on the synchronization dynamics is more pro¬ 
nounced than that of body translation [inset, curves (hi) 
and (iv)]. Without body rotation, (9) converges more 
slowly towards zero, qualitatively resembling the results 
of anchored flagella [26]. However, absence of body trans¬ 
lation does not affect the outcome much. Hydrodynamic 
interactions (HI) between flagella are important in repro¬ 
ducing the standard result of the full simulation [inset, 
curves (i) and (ii)]. However, we find that synchroniza¬ 
tion is even possible without HI contrasting the situation 
of anchored flagella, where HI between flagella are known 
to be essential [26]. Interestingly, cell-body movement is 
sufficient to synchronize flagella, similar to findings for 



t [ms] 

FIG. 4. Comparison of different factors affecting bundling 
and synchronization dynamics (inset). Curves (i): standard 
simulation run performed for = 23°, 9{i) = 0.8 Vi, and 
Tm = 3.4pN/rm. Other curves show simulation runs that 
ignore: (ii) hydrodynamic interactions among flagella, (iii) 
axial body rotation, (iv) body translation, (v) steric interac¬ 
tion between flagella, or (vi) both steric and hydrodynamic 
interactions. 


Chlamydomonas in Ref. [39]. 

Bundling dynamics is affected with a similar trend. 
While bundling is delayed by ca. 10 ms in absence of HI 
between flagella, absence of body translation is less se¬ 
vere [curves (i), (ii), and (iv) in Fig. 4]. There is no 
bundle formation in absence of axial body rotation [ hor¬ 
izontal line as curve (iii)] because the body re-orients and 
slows down translation, resulting in an unusual buckling 
of flagella away from each other. More significantly, when 
we follow Ref. [30] and allow temporal evolution of each 
flagella affected only by body movement but not by ei¬ 
ther HI or steric interactions, bundling gets significantly 
delayed by about 30ms [curves (v) and (vi)]. 

To conclude, we present a detailed modeling of an E. 
coli with two flagella and a motile cell body. This al¬ 
lows us to probe bacterial propulsion on an yet experi¬ 
mentally inaccessible level. The complex role of flagellar 
polymorphism is ignored for simplicity. In principle, it 
can be probed extending our model ]17]. We demonstrate 
that compared to the situation of anchored flagella, flag¬ 
ellar dynamics close to real conditions is strikingly altered 
by body movements. Times to bundle and synchronize 
are dramatically reduced. In particular, flagellar entan¬ 
glement helps bundling to proceed quickly like a ‘zip¬ 
ping’ motion, which we rationalize in a simplihed model. 
Furthermore, we demonstrate that body movement and 
flagellar entanglement lead to rapid bundling and syn¬ 
chronization even when hydrodynamic interactions are 
neglected. Our findings are important in explaining ex¬ 
perimentally observed times scales as mentioned above. 

Finally, more and more artificial microswimmers using 
different swimming mechanisms have been and are con¬ 
structed [40] . We provide here an example how one devel¬ 
ops a model for exploring and ultimately understanding 
the biomechanics of microswimmers. 
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